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ABSTRACT 

The extended wide-angle parabolic wave equation applied to electromagnetic wave 
propagation in random media is considered. A general operator equation is derived 
which gives the statistical moments of an electric field of a propagating wave. This 
expression is used to obtain the first and second order moments of the wave field and 
solutions are found that transcend those which incorporate the full paraxial 
approximation at the outset. Although these equations can be applied to any propagation 
scenario that satisfies the conditions of application of the extended parabolic wave 
equation, the example of propagation through atmospheric turbulence is used. It is 
shown that in the case of atmospheric wave propagation and under the Markov 
approximation (i.e., the (5-correlation of the fluctuations in the direction of propagation), 
the usual parabolic equation in the paraxial approximation is accurate even at millimeter 
wavelengths. The comprehensive operator solution also allows one to obtain expressions 
for the longitudinal (generalized) second order moment. This is also considered and the 
solution for the atmospheric case is obtained and discussed. The methodology developed 
here can be applied to any qualifying situation involving random propagation through 
turbid or plasma environments that can be represented by a spectral density of 
permittivity fluctuations. 
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I. INTRODUCTION 


Since its inception by Leontovitch and Fock 1 , the parabolic wave equation has 
found application in all aspects of electromagnetic wave propagation scenarios 2-4 as well 
as those in acoustic wave propagation 5 " 7 . In the case of electromagnetic wave 
propagation in random media, the parabolic wave equation was originally applied by 
Klyatskin and Tatarskii 89 to describe the statistical moments of the attendant electric 
field. This type of parabolic wave equation, however, is accurate only to within the 
paraxial approximation and is therefore suited to treat small-angle wave scattering about 
a preferential direction of propagation. Thus, in the case of propagation through random 
media, the parabolic wave equation can strictly be applied only when the smallest size / 0 
of the permittivity fluctuations within the medium and the wavelength A satisfy the 
condition /„ » A, i.e., the wavelength must be the smallest spatial scale in the problem. 
In the case of the open atmosphere, ^-Imm. As shown by Klyatskin and Tatarskii 10 , this 
is indeed a good approximation in the case of atmospheric propagation at optical 
wavelengths. However, there are situations in the applications, e.g., atmospheric 
millimeter wave propagation, in which one has / 0 a A and the paraxial approximation 
strictly no longer holds. In these cases, one can have wide-angle scattering about the 
preferential direction of propagation. Again, in the scenario of electromagnetic wave 
propagation through random media, such considerations have been made 11 . Much more 
recently, the same situation has been met withyn acoustical wave propagatioi^j 1 ^ 15 . In 
what is to follow, the wide-angle modification of the parabolic wave equation will be 
called the extended parabolic wave equation for reasons that will be apparent in the next 
section. Although extensions of this theory have been advanced within the realm of 
stochastic propagation through random media 11 ' 16 , none have been analytically treated at 
the level where corrections to the paraxial approximation have been quantitatively 
identified and compared to paraxial results. 

It is the purpose of this work to analytically derive solutions for the first and 
second order moments of a plane wave field propagating through a random medium 
within the context of the extended parabolic wave equation. Approximations must, of 
course, be made in order to obtain analytical results. However, the main point here is that 
such approximations are not made at the outset, as is usually done. Such solutions will 
afford a comparison between the results of the parabolic wave equation within the 
paraxial approximation and those of the extended wide-angle theory. In Section 2, the 
extended parabolic wave equation for electromagnetic wave propagation is derived, for 
completeness, from the stochastic Helmholtz equation. Since this will be an operator 
equation in the random electric field, a statistical operator method is developed in Section 
3, which will give a general equation for an arbitrary spatial statistical moment of the 
wave field. From this, expressions are obtained as special cases for the first and second 
spatial moments of the field. The first field moment is just the average field and the 
second field moment is known as the mutual coherence function (MCF). Analytical 
solutions to these equations are derived for the Kolmogorov spectrum of atmospheric 
permittivity fluctuations within the Markov approximation. Unlike earlier treatments of 
propagation using the parabolic wave equation 8 " 10 , the generality of the operator method 
also allows one to obtain an expression for the second order field moment in the direction 






longitudinal to the direction of propagation; this is known as the generalized MCF. A 
solution for this quantity is also obtained. 


II. THE EXTENDED PARABOLIC WAVE EQUATION FOR RANDOM MEDIA 

Consider the scalar stochastic Helmholtz equation for an electric field E(x,p) 
propagating principally along an x-axis and perpendicular to the p -plane of an otherwise 
arbitrary coordinate system 


d 2 E(x,p) 

dx 2 


+ V 2 £(x,p) + k 2 E(x,p) = k 2 e(x,p)E(x y p) 


( 1 ) 


where e(x,p) is the random part of the total permittivity e(;c, p) = 1 + e(x,p) of the 
propagation medium. Using the decomposition of the total field into a forward 
propagating field E + (x,p) and a backward propagating field £~(x,p), i.e.. 


E(x,p) « E + (x,p) + E~(x,p), 


3E{x,p) _ dE*(x,p) dE~(x,p) 
dx dx dx 


(2) 


with the expansion into inhomogeneous plane waves 

£*(x,p)=jf / e*(q)expiq\p±i(k 2 -q 2 f 2 x\d 2 q 

—GO 

one can obtain the following operator expressions for the component fields 


(3) 


JE ± (x,p) 


2i ^^ ±1 ^ +K f E , [xrp)= ^ k , + ^ 


(4) 


Such an equation (or its Fourier transform) was considered by Malakov and Saichev 17 as 
well as by Klyatskin 18 and Frankenthal and Beran 19 . 

In the event that the backscattered field is insignificant with respect to the forward 
scattered field, i.e., E + (x,p) »E~(x, p) , one can formally set E - (;t,p) = 0 ; Eq.(4) thus 
becomes a single equation for the forward propagating field 

This equation, which can be called the “extended parabolic equation ,,)S , was previously 
analyzed by Saichev 16 using a method different from the operator formulation in w T hat is 
to follow. Its integral formulation is also known as the method of multiple forward 
scatter 11,20 . However, it was pointed out 11 that this method is applicable to situations 



when 21 A £ / 0 , the penalty for ignoring the backscattered component; The classical 
parabolic wave equation in the paraxial approximation holds for cases where A << l 0 . 
Hence, although Eq.(5) is capable of describing propagation situations in which the wave 
is scattered at angles up to Jtj 2 with respect to the x axis, its application is limited to 
cases where A ^ / 0 . 

The form of Eq.(5) can be simplified by defining the differential operator 




( 



1/2 


( 6 ) 


and the corresponding integral operator 


Ik 2 


l v*\-^ 

1+ / 

\ K , 


Equation (4) then becomes 


2ik '^x^‘ + U p E ( x 'P) - 2kAy p*( x ,P) E { x >p) » °, 



(8) 


(Hereafter, the superscript ‘+ will be dropped.) This equation can be reduced to the 
well-known stochastic parabolic wave equation in the paraxial approximation by 
expanding the operators to give U p ~2 k 2 + V‘ and V p ^ijlk 1 and transforming the field 
via E(x,p ) = W(x,p)exp(/£x) . It must be noted that these operator expansions, which 
are necessary conditions for the paraxial approximation, applied to Eq.(8) at this juncture 
yields the paraxial form of the parabolic wave equation. In what is to follow, the full 
form of Eq.(8) will be dealt with in forming an expression for the statistical moments of 
the electric field. Only during the solutions of these equations for the first and second 
moments will these operator expansions be judiciously applied. It will turn out that a 
geometrical optics approximation which is implicitly made in the scattering term of these 
solutions can be lifted which will essentially allow diffractive effects to be introduced 
into the scattering term. 

In what is to follow, Eq.(8) will be employed to yield an operator equation for the 
generalized statistical moments of the field. Here, the term ‘generalized’ connotes 
moments at differing transverse coordinates p as well as differing longitudinal 
coordinates x . ■ . 


III. OPERATOR SOLUTIONS FOR THE GENERALIZED FIELD MOMENTS 
A. An expression for the generalized nm* field moments 

Defining the stochastic operator 
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D ,.p- 2ik -^+(J p + zk \K x ’P) (9) 

Eq.(B) simply becomes 

D ,. p E{x,p)=0 (10) 

which is easily amenable to further statistical analysis. To this end, one defines the 
generalized moment of the electric field 

^r™{ X i'Pl’ X 2’P2’ m '’ X n ’Pn’ X n+rP n +v'' X , l+m iPn+m) -{gnm}' (H) 


n+m 

8~ m ]j E { x j'Pj)n E *{ x i*pi)- 


( 12 ) 


j - 1 


f-n+l 


Using a modification of previous prescriptions 22 23 , one can employ Eq.(10) for each of 
the field points [x, p) and obtain for the product of fields 


L n m8 n m( X l’A’ X 2>p2'‘‘‘ X n’Pn> X n+vPn+i'' m ' X n +m'Pn + rjt) = Q 


where 


L = 2 ik 

run 


i n+m ~ \ n 

a x-\ a 


{ 2^2^2K^M)- 

\;-l l-n+l OX l ] y -1 


!-n + 1 


In order to isolate the quantity (g wn ) = T^from this relation, it is expedient 24 ' 25 to 
decompose the operator L nm into its average and random parts, i.e.. 


(13) 


n+m 

-2 K+2*%f(x,. Pl )) (14) 


n+m o 


U-1*I dX i J M 


/-n+- 1 


(15) 


and 


- 2 k‘ 


n n+m 

2 V (• x j’Pj)~ 2 K r ( x „p,) 


i-\ 


/“«+i 


’ ( Z -) = ° 


06) 


Hence, Eq.(13) becomes 


((O + 4*,)&w =0 * 


(17) 



Ensemble averaging this relation yields 

{ Aim + (4ra^nm) “ 

Similarly writing 

Snm ~ ^ Snm’ { Snm ) — 

and substituting into Eq.(18) gives 

(^)r |w , + (j-'iun8wn) ~ 


(18) 


(19) 


(20) 


Remembering that it is the goal of this development to obtain an expression for the 
general field moment T nm , one subtracts Eq.(20) from Eq.(13) and using Eq.(19) obtains 25 

+ ^nmSnm ~ *“ ran Snm) = 0 . (21) 

Combining the first and third members of this equation using the fact that 
-(0] r - - gives 

+ 0. (22) 


One must now isolate the random quantity i g nOT by defining an operator L n l m inverse to 
>•©•» L ~,l L nm = 1- Thus, operating on Eq.(22) with L x nm yields 

l^ + ^i^r wit -Li(L nm | nm ) = 0. (23) 

Finally, operating on this relation with L nm> ensemble averaging and solving the resulting 
expression for gives 


{^rnnSnm) [l • (24) 

Substituting this result back into Eq.(20), one obtains for the equation governing r nCT 

{(^>-[i-(^C)f(£™C,£™}}r m ,=o. (25) 

The solution of this operator equation gives an exact solution for the arbitrary field 
moments for wide-angle propagation through a random medium characterized by the 
stochastic permittivity e(x,p) and the assumption that ^£^ = 0. 



The general relation given by Eq.(25) can be reduced to the parabolic equation for 
the field moments in the paraxial approximation in the case where A «/ 0 . In this 
instance, one employs the approximations for the operators U p **2 k 2 + V* and 
V p “ l/ 2A: 2 used earlier. In addition, the classical parabolic equation considers statistical 
moments in the same transverse plane, i.e., x, ~x, = Thus, the partial differential 

operators in Eq.(15) collapse into the single operator d/ax. Equations (15) and (16) then 
become 


n n + m 

^> 2ik i+lK- iv 


(26) 


H 


I-l+n 


and 


Km 


« n+m 


i J-l 


l-l *n 


(27) 


Two related palatable approximations must now be made; since s-(l \ + L . and it 

ntn \ nw/ ww ’ 

is usually assumed that ||Z nm [|«l, one has 


where 


- [{L m ) + L m f - (£ J", [l - (L m L^ )]■' 


(28) 


(U"' - 


jn n n+tn 

2S-+\v- -Yv 
& p ‘ A 


j-l l-l+n 


1-1 


2 ik— 
dx 


1-1 


« — ]dx\ 
2 ik i 


(29) 


Equation (25) then becomes 


2 ik~+yv 2 p -Yv* -lL (L ] )L ) 

fc Zj pj jLi P> \ «» \ W" / nm / 


/-l J-l+n 


1^=0, 


(30) 


where 


22 fa’pMxp/))- 2 


v-if-i 


j-l f-1 +n 
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rt (Mm ' n4 m n+ m 

~E 2 [z( x ’PM x '-P'>)) + E E (e*(jc,pi)s*(jc',p;.)) 


;’-l M+« 


i-i+fl r-i+B 


<*4 (31) 


which is the well known classical paraxial form for the problem 23 . It is interesting to note 
the ‘geometrical optics’ approximation made in Eq.(29) which gives the usual scattering 
term of £q.(31). One can envision a substantial extension of this development beyond 
that of the classical treatment if one is to use the entire form of the operator , i.e., 

use the inverse of the operator ( L mn ) as solved in its entirety rather than in the 
geometrical optics approximation which was used above. The inclusion of the V 3 
operators will allow diffraction phenomena to be incorporated into the scattering term. 
This will form the subject of the next section in which the first order moment and the 
generalized second order moment (i.e., the generalized MCF) are derived using the above 
formalism. 


B. Solution for the first moment 

The first order moment or average field of the random electric field in a plane transverse 


to the direction of propagation is defined through Eq.(ll) to be given by 

r .o(*;A) = (£(*;A)) (32) 

which is a solution of the operator relation of Eq.(25), viz, 

{<Ao>-[i-(i,oAo)]" ! (A^„)}r 10 = o (33) 

where, from Eqs.( 15) and (16), 

{i^o) = 2ik-?- + U p , Z, 0 =2k A V p l(x,p), (34) 

The ability to proceed in an analytical fashion is dependent upon the simplifying 
approximations used earlier, viz,, one has that 

^io ~ [{^io) + Ao] ^(Ao) » [l-(AA"i)] ' as i* (35) 

Using this result, Eq.(33) becomes 

|{Ao) - ^Ao{Ao) £o)}r,o- 0 . (36) 

Employing the appropriate definitions of the operators, Eq.(36) gives 


R 



-1 


m K + u * '((^ ^ (u 4 vX*'P)) 


r io (^p) = 0. (37) 


This differential equation in the operators U p and V p must now be simplified and solved 
for the first-order moment F 10 (x, p) . 

To this end, one must first deal with the ‘diffraction propagator’ 


2 ik—+U p 
dx p 


25 G{x,p) 


(38) 


It is at this point that this method has its greatest effect on the description of the 
propagation mechanism. Within the geometrical optics approximation, the U p term in 

Eq.(38) is neglected with respect to the Hkd/dx term. Hence, potentially important 
diffractive effects are neglected within the scattering term (• •) of Eq.(37) by use of this 
approximation (see Eq.(29)). As mentioned earlier, the retention of the U p term, even in 
its approximate form, will include diffractive phenomena inherent in the random 
scattering mechanism that have hitherto been neglected. The Green function of the 
operators 2 ik d/dx + U p given by Eq.(38) is defined by 

G(x,p) = d(x - x')6(p -p') (39) 


d 

v 2 ;l 

V- 

2ik~ + 2k 2 



dx 

k 2 

K ) 



where the definition of U p is used. Applying the approximation U p «2 k 2 +V 2 and 
solving for the Green function G(x,p) yields. 


r e X p[-i*(p-p') 72 (*-* r ) 

p[-z£(x -x)\ ; . 

x-x 


G(x,p) = G(x,p;x\p') = [ — ex 

\4jtJ 


Thus, the third term within the brackets of Eq.(37) can be written 
{->■ e{x,p))\^kj-+ t/,j ‘(2 kXK*'p)ij = 

X 00 
0 -® 

Proceeding further, one now must deal with the operator products 


(40) 


(41) 
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V^pj-j 


( 1 \ 

V 2 


l + ~f 

\2* 2 j 

k 2 

K 7 


-1/2 


r(x,p). 


(42) 


Since s{x,p) is a random function, it can be represented in the form of a Fourier-Stieltj 
integral 36 , i.e., 


es 


e(x y p) = f exp {iK-p]dZ(x,K) 


(43) 


in which the spectral amplitude dZ(x,K ) is endow' ed with the same statistical properties 
as is the random function e(x, p) as will be shown in what is to follow. Applying Eq.(43) 
to Eq.(42) gives 


V p e(x,p) = 


2 k‘ 


I/-F 




exp(/k* p)dZ(x,K). 


m 


Thus, the ensemble averaged product appearing in right side of Eq.(41) becomes. 


{V p e(xrp)V/e(x\p'))J^\ ff 




V ^/2 


i-^r 

k 2 


•exp (iK-p + irc''p')(dZ(x,K)dZ(x\fc')). (45) 

One now makes use of the fact that the atmospheric permittivity fluctuation field s(x,p) 
is taken to be statistically homogeneous, characterized by a power spectral density <J> e (jjf) 

and d-correlated in the longitudinal direction (the Markov approximation); these 
circumstances allow one to write 26 


(dZ(x ) fc)dZ(x , ,K , )) = d(/c+ k')F s (x - x\K)d 2 Kd 2 K' (46a) 

where for d-correlated fluctuations in the x direction, the two-dimensional spectrum 
is given by 

F e (x - x\ k) = 2jtd{x - x')$> E (k) (46b) 

m which is the three dimensional spectrum of permittivity fluctuations. Using 

these relations in Eq.(45) and performing the integrations where possible yields 

= d(x-or')/ exp [iK-p d ]<P E (K)d 2 K 

(47) 


10 



where p d ~ p-p is the difference coordinate. 

Taking the statistics governing the random field e(x, p) to be also isotropic, 
i.e., <£,(£) Eq.(41) can now finally be evaluated by substituting into it Eqs.(40) 

and (47); converting the integration in the p^-plane into one in plane polar coordinates 
and performing the associated integrations gives 



~ i3tV f [ J- fr 


-1 

<P £ (K)KdK 


(48) 


where the d-function relation 

fd(x-x')dx '-t 
0 z 

is employed. 

Returning to Eq.(37) and, substituting Eq.(48) into Eq.(37), one obtains 


(49) 


2ik~+2k 2 

dx 


( 

K / 


L-I 


+ at 2 A - 3 f l-^j d> K ( K ) K dK 

n l k , 


r 10 (^p)=o (50) 


The general solution to this equation is unknown. However, for the plane- wave case, one 
has that 


( vV 

V2 

f v 2 \ 


r io(*.?) = 

1+ # 


W - r,oW 


(51) 


since the plane wave will not possess any transverse variations. In this special case, 
Eq.(50) becomes 


2ik~ + 2k 2 + ijr*k 3 f lUp-j *'{K)KdK |r io (x) = 0 


dx 


(52) 


the solution of which is 

r ioW - r io(°)exp 


, 2 cc 


>-1 


ikx-~k 2 x§ 1 - <3> e [K)KdK 

2 n k 


(53) 


This result differs by the factor ^1 ~ k 2 j , k 2 ^j from that of the parabolic equation in 

the paraxial approximation. Of course, the later is obtained from the former by retaining 
the first term in the series expansion of this factor. The presence of this factor tends to ° 



accentuate the spatial frequencies near the value of the wave number k. Since this theory 
is applicable to those situations in which As 7 0 , most of the contribution of this factor to 
wave scattering will occur at the largest spatial frequencies of the inhomogeneities. 
Hence, in the case of atmospheric turbulence, one can employ the von Karman spectral 
density (i.e., the Kolmogorov spectral density modified by a lower cuf off frequency K 0 ) 

®t( K ) = 0.033C c 2 (k 2 + (54) 

which is not bounded at the high spatial frequencies, to compare the result of Eq.(53) 
with that of the paraxial approximation. Here, C] is the ‘structure parameter’ giving the 
strength of permittivity fluctuations and K 0 = 2 JtjL 0 where ^ is the outer scale (or the 
largest spatial size) of permittivity inhomogeneity. Substituting Eq.(54) into Eq.(53) and 
evaluating the integral using Mathematic a 21 yields 


r»( ;i: ) ~r io (0)exp 


ikx 


(1.772+1.023/)*-^ + 0.3908A' 0 ^ 2 /^ 


’V k 2 


>C:k 2 x 


(55) 

where the structure parameters for the permittivity C] and the associated refractive index 
C n are related by C 2 = 4 C 2 . In the case of the open atmosphere, one always has within 
the bounds of the extended parabolic wave equation, k»K 0 ; thus, the hypergeometric 
function reduces to unity and k^ 3 « Kf 3 , allowing Eq .(55) to be approximated by 


r,oM - r io (0)exp[«x -0.3908A'- 5 ' 3 C a 2 £ 2 *] (56) 


which is the result of the parabolic equation in the paraxial approximation. Hence, the 
use of the extended parabolic wave equation only makes negligible amplitude and phase 
corrections to the first order moment of the wave field propagating through atmospheric 
turbulence. This result extends and establishes the accuracy of the parabolic equation in 
the paraxial approximation for the first order moment (mean field) as it applies to 
atmospheric turbulence at all acceptable wavelengths of application, A-s / 0 . The next 
section will consider the calculation of the generalized second-order moment (generalized 
MCF) of the wave field from the extended parabolic equation and compare its result to 
that of the paraxial approximation. In addition, due to the completeness of the operator 
analysis, one naturally obtains expressions for the MCF along the longitudinal axis. 


C. Solution for the generalized second moment 

The generalized second order moment or generalized MCF of the random electric field in 
two planes transverse to the direction of propagation is, from Eq.(l 1), 

^"j]('*'l»Pp , *2’P2) = {^(-^1’Pl (*^2’P2)} 


1 ? 


(57) 



which is a solution of Eq.(25), in this case given by 

{(A .) - f 1 - (A )] (A ,)}r. ■ - ° 

where, from the definitions of Eqs.(15) and (16), 


(58) 


and 


(Ai) = 2/i— + 2ik—+ JJ -V' 
' ' ax, ax, 


^ s2ki \yA x ^)-V‘ p e{x 2 ,p 2 )], (A,} = 0. 


(59) 


(60) 


As with the case for the first order moment, two related approximations must be made at 
the outset to render the problem analytically tractable. In particular, so long as III, J «1, 


A!-[(Ai) + Ai] ”{Ai)\ [i-<AiA!)] 

which allows Eq.(58) to become 

{(Ai)“(Ai(Ai)A,)}a 1 =0; 


(61) 


(62) 


At this point, it is suggested to connect the longitudinal coordinates x, and x., to the 
related centroid x c and difference x d coordinates. 


x, + x~ 

x 4 mx l -x 2 . 


(63) 


The operator expressions of Eqs.(59) and (60) then become 




A, -2 + -f- .ft)]. <4> = o. 

Hence, Eq.(62) can be written as 


(64) 


(65) 


n 



2 ik—+U-Ul -l\2k 4 \v p {x c + ^.p}-V’e'(x-^-,p 2 

dx. p ' P! \ M e 2 p ‘ V 2 2 


1 * + °* VK*' t’A W' K- ) r n -°- 


The solution of this equation commences with obtaining an expression for the 
Green function, i.e., the diffraction factor 


2il * * U »-K>\ =G{x c ,p v h). 


Proceeding as in the last section and using the approximation U p «2 k 2 +V 2 p , this 
requires the solution of 


m ^T + V -' V - | G ^A>P 2 ) -x’Mpi - pMPi -Pz) 


which is given by 


ik\( i w i 




= G(x c ~x' c ,p l -p[,p 2 -p 2 ) 


Therefore, the fourth term in Eq.(66) becomes 


(•••> = (2* 4 ) - / f J G(x,p,,p 2 ;x'p;,p;)^y Ci v;.E^ c + ^,p 1 js^>^,p 1 'j- 


- v „ v k i *c + V’P 4-^-.Pi -V^g* ir^pJt ^fip; + 


+v lylf (vo) 


One now employs the Fourier-Stieltjes transform as before to represent the products of 


V ’ Le -> 


14 



X , - 


V *c * TT»A. 2 = 


a 

2fc 2 


/ 




-yi 


exp(/K, 2 • p i 2 )c/zl x c ± ^- s ic 12 j, etc., (71) 


and obtains the following relations (functional arguments have been suppressed but the 
correspondence to those in Eq.(70) follows) 


(w g )- 


1 


\2 


-2V 1 


2k‘ 


f 1-^- expfiif, ■ (p, - pOKF - x 'cA) d \ * (72a) 


/(('fr) exp [ i ^-{P,-Pi)]F,{ x c-K + x J A) d \- ( 72b ) 




2 ( K 2 ^ 

k 2 


\ 


exp[-i k 2 • ( p 2 - A')] ( x c ~x c -x d , K 2 )d\ , (72c) 


(«">- 3F S 


1 K* S 
\ */ 


exp[-iic 2 • (p 2 -#)]F,(* e ” x;,ff 2 )j 2 /c 2 . (72d) 


Substituting Eqs.(69) and (72a)-(72d) into Eq.(70), performing the required integrals over 
A and p 3 in plane polar coordinates, and taking all the spectra /^(x.?, 2 ) = F f (x,fC l2 ), i.e., 
to be isotropic in the frequency k, one obtains 




0 0 




.Ax c -K'i 

N 

exp 

l ' ^ ^ • 

2k 


-exp 

-exp 


y&'-K) 

2k 

-Xc) 


-l 


2k 

+exp 


J 0 ( k p*) F '( x c- x ' +x <’ k )- 

■lo{ K P^( X c- X c- X i’ K ) + 

KdKdx„ 


.K 2 (x c -x' c ) 

F f {x. -x’k) 

2k 

E \ C C 7 / 


(73) 


where p d s|p 2 - p,|- Using Eq.(46b) in Eq.(73) and performing the x c integration, 
remembering Eq.(49), finally gives 


{•••) = -2 x 2 ik 3 J 

0 

Hence, Eq.(66) becomes 


f ) ,2\' 1 

1- Fi 


1 - exp - 


■ 2 \ 
IK X ' 


2k 


Jo{ K Pj)\®e{ K ) KdK - ( 74 ) 



2N 


2 ik 


Ik — + U p - U' +2x 2 ik 3 f 

* P2 'f. * 2 


-i 


1-exp - 

2k 


j o{kp<) 


4> f (/f 


f„-o 


( 75 ) 


where r, } - T, ,(*c»*,mP*)* Making the plane wave approximation of Eq.(75), analogous 
to that done earlier for Eq.(50), the resulting differential equation has as a solution 


r(x c » , p d ) — r (0, 0, p d ) exp 


-x 2 k 2 x,.f 




i-£ 


o \ 


1- exp - 


iK 2 x. 


2k 


Ko(«Prf); 


-® s (k)kcIk 

(76) 


It is important to note the initial condition r(0,0,p rf ); since one necessarily must take 
x c =0, one must also have - 0 by Eqs.(63) since jc,,jc 2 & 0. 

This result cannot be analytically studied in its entirety and hence will be 
considered in two special cases. The first case is defined by x d =0 in which one deals 
with the usual transverse MCF T(x c ,p d ) = r(x c ,0,p d ). Thus, the extended parabolic 
equation solution for the MCF in a transverse plane at a distance x c from the source is, 
fromEq.(76) 


r i x c’Pt ) = r(0,0, p„)exp 


-x 2 k-xj 


l-Sr 


0 \ 


•) )}«>«(*•)*“** 


(77) 


As with the case of the first-order moment, this result differs from that of the paraxial 
approximation in the presence of the factor (l-* 2 /* 2 )" 1 which serves to accentuate 

spectral contributions for frequencies near k. The integral indicated in this expression 
cannot be analytically evaluated using the von Karman spectrum, Eq.(54). However, 
since the integrand is such that no singularities exist in the use of the unbounded 
Kolmogorov spectrum, viz., Eq.(54) with K 0 =0, one can use such a spectrum in Eq.(77) 
and, upon evaluating the integral via Mathematical obtain 


r n(-V>Prf)^ r n(°’P</)expi 


-j(3.544 - 2.048/)r^(l-/ o (^ t/ )) + 


(78) 


where 1 F 2 (---) ’ s a generalized hypergeometric function and, as noted earlier, C 2 = 4C 2 . 
At the outset, since A « p d in most applications, the first term within the braces of 
Eq.(78), although quite interesting in structure, is negligible with respect to the second 
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term. The hypergeometric function of the second term is most easily dealt with by first 
converting it to a Lommel function 28 .*„<,(•••), i.e.. 


1^2 


Lij 17. k*p) 
’ 6 ’ 6 ’ 




yj (wr% 0 fe) 


(79) 


The Lommel function reduces 28 , in the case where kp d » 1, to the simple approximate 
result sy 0 {kp d )~(kp d f\ Using this in Eq.(79), Eq.(78) becomes the well-known 
paraxial result 


^u{ x c^Pd)^ r u(^Pd)^v[-^^^Clpfx^ 


(80) 


Again, as with the first-order moment, corrections to the second moment afforded by the 
extended parabolic wave equation are negligible in the case of atmospheric turbulence. 

The second special case in which Eq.(76) will be examined is in the instance 

where the factor (l-tf 2 //: 2 ) can be neglected; as shown above, this is a good 

approximation for atmospheric turbulence. Thus, one is now dealing with the solution 
for the generalized MCF for field locations at different transverse and longitudinal points, 


r(* c » x d ,p d ) = r(o,o ,p d )exJ -n 2 k\] 


1 - exp - 


iK 2 x d 


2 k j 




. (81) 


Longitudinal correlations of the wave field have been previously considered 29,30 using 
different methods and obtaining different results. Setting x d = 0 in Eq.(81) gives the well 
known paraxial result for the transverse MCF r(x e ,0,p d ). Unlike the cases studied 
above, the form of Eq.(81) suggests the use of the Kolmogorov spectrum as modified by 
Tatarskii, which incorporates a cutoff at high spatial frequencies by allowing the 
introduction of the inner scale of turbulence / 0 , viz.. 


«>,(*) -0.033C„V ,,3 exp _ 
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2\ 


K, 
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m 1 


5.92 
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(82) 


Substituting Eq.(82) into Eq.(81) and evaluating the resulting integral yields for the 
generalized MCF 


rfo, x d , p d ) - r (0, 0 , p d ) exp 


-A352k 2 Clx c \ B 3,0 ,F,| 




Pd 


6 4 B 
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(83) 


where 
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_ 1 ,X A 

B ^ + ‘f k m 

Thus, as previously noted 30 , the presence of the diffraction factor on Eq.(83) modifies the 

effect of the cutoff frequency K m . In the case where p d >\rfjc 2 m + ix a /2f^ 2 , Eq.(83) 

reduces to Eq. (80) upon employing the asymptotic representation of the confluent 
hypergeometric function ,/*[(*••). When p d =0, Eq. (83) becomes 



^,o)=r(aap d )exp 
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. 2 ^ 
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(85) 


In the appropriate limits, this expression gives 

f exp[-3 .8 1 33 i kC 2 Kfx c x d ] , 


rU c ,*,,o) = r(o,o,p d )■ 


x d k«ll 


exp[-(0.563 1 + 2. 1 02 z )k lf6 C 2 n x c xf ] , x d k » / 0 2 


(86) 


Hence, there is a phase variation in the longitudinal direction, as expected, with an 
attendant attenuation as the longitudinal separation is increased. It must be noted, 
however, that one needs to realize the condition x d > l 0 in all cases so as to satisfy the 
assumption of d-correlation of the fluctuations along the longitudinal axis 30 . 

IV. SUMMARY AND EXTENSIONS OF THE FOREGOING 

A general operator equation for the moments of an electromagnetic wave field 
propagating through a random medium as described by the extended parabolic wave 
equation, i.e., Eq.(5), has been derived and is given by Eq.(25). Unlike the usual 
parabolic equation approximation in the paraxial approximation which holds only for 
A << / 0 , the extended version of this equation holds for cases up to As / 0 ; this limitation 
is dictated by the neglect of the backscattered field since the extended parabolic wave 
equation describes scattering at angles up to jr/2 about the preferred direction of 
propagation. An example application of this theory was chosen to be that of the scenario 
of propagation through atmospheric turbulence in which the first and generalized second 
order statistical moments are given by Eqs.(50) and (75), respectively. In the course of 
deriving these equations, use of the paraxial approximation is avoided at the outset, 
allowing one to be able to include diffractive effects within the scattering terms of these 
equations. These relations, which assume a <5-function of permittivity fluctuations in the 
preferred direction of propagation (i.e., the Markov approximation) are then solved in the 
plane wave case yielding Eqs.(53) and (76) in terms of the spectrum & e (k) of 
homogeneous and isotropic fluctuations. Kolmogorov type spectra (i.e., spectra ~ /E 1 * 3 ) 
are employed to obtain analytical expressions for the first and generalized second 
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moments. In this atmospheric example, this extended theory showed that the paraxial 
approximation holds even for millimeter wavelengths in which A ~ / 0 . This is due to the 
fact that the spectra which is used to represent atmospheric turbulence (all variations of 
the Kolmogorov spectrum for atmospheric turbulent fluctuations) is such that the 
contribution of the spectral frequencies, which approximately correspond to millimeter 
wavelengths, is relatively small as compared to those at the smaller spatial frequencies 
(which correspond the outer scale of turbulence). This gives rise to the insignificant 
levels of the obtained correction terms at nominal operating wavelengths. However, the 
novel feature of this application example is the description of the longitudinal second 
moment which naturally follows from the comprehensive operator solution and the 
incorporation of diffraction within the scattering term. A solution for this quantity, given 
by Eq.(83), is discussed as well as its limiting forms. 

The operator method presented here from which the statistical moments are 
obtained is general enough to allow the use of assumptions less restrictive than the 
Markov approximation. The use of fluctuation spectra that do not decay as rapidly as that 
of the Kolmogorov case (e.g., that of a turbid medium or some turbulent plasmas) will 
give rise to non-negligjble correction terms for the first and second moments of the wave 
field. 
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